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FORWARD 

This  report  describes  work  done  at  the  Gas  Turbine  and  Plasma 
Dynamics  Laboratory  (GT&PDL)  at  the  Massachusetts  Institute  of 
Technology  for  the  Air  Force  under  Grant  No.  AFOSR-83-0034 .  Dr. 
Anthony  K.  Amos  was  the  technical  monitor. 

The  work  reported  herein  was  performed  during  the  period  1 
December  1983  through  31  December  1984.  Three  graduate  students,  Otto 
W.K.  Lee,  Ya-Pei  Chang  and  Eva  Czajkowski  worked  under  the  supervision 
of  Professors  Manuel  Martinez-Sanchez  and  John  Dugundji  at  some  time 
during  this  period. 
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1.  INTRODUCTION 

As  explained  in  the  previous  Yearly  Report,  the  purpose  of  this 
program  was  to  extend  our  knowledge  of  the  fluid  forces  arising  from 
the  presence  of  labyrinth  seals  in  high  power  turbomachinery,  and  to 
relate  these  forces  to  their  dynamical  effects  on  the  turbomachine 
structure.  Thus,  the  investigation  was  divided  into  two  parts,  one 
dealing  with  the  fluid  aspects  and  the  other  with  the  structural 
dynamics  part  of  the  problem.  This  report  will  detail  our  progress  in 
each  of  these  two  areas  during  the  second  year  of  our  Grant. 


2.  WORK  PERFORMED 

2.1  Fluid  Dynamics  of  Seals. 

(a)  Test  Facility 

In  previous  the  Yearly  Report  (Ref.  1)  we  described  the 
formulation  of  an  analytical  model  of  the  fluid  mechanics  of  a 
multi-chamber  labyrinth  which  operates  in  the  presence  of  shaft 
whirling  motion.  A  preprint  which  was  later  presented  (Ref.  2)  at  the 
3rd  Workshop  on  Rotor  Dynamic  Instability  Problems  in  High  Performance 
Turbomachinery  (Texas  A&M  University,  May  28-30,  1984),  was  included  as 
an  Appendix. 

One  of  the  important  points  that  came  out  of  this  work  was  the 
fact  that  theory  predicts  seal  whirl-dependent  (damping)  forces  of  the 
same  order  as  the  static  (stiffness)  forces,  and  yet  no  experimental 
data  were  available  for  validation.  This  has  led  us  to  undertake  the 
design  of  a  seal  test  rig  in  which  the  rotor  will  have  both  spin  and 
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whirl  notion,  and  these  will  be  controllable.  Since  theory  also 
indicates  that ^non-pressure  forces  are  secondary,  measurements  will  be 
made  of  the  time-dependent  non-uniform  pressure  pattern,  to  yield 
cross-forces  by  integration. 

In  addition  to  the  ability  to  produce  controllable  spin  and  whirl, 
the  requirements  on  this  test  rig  include  (a)  Realistic  pressure  drops, 
i.e.,  a  pressure  ratio  of  1.3  in  the  last  seal  strip.  (b) 

Constriction  Reynolds  no.  of  at  least  1-2x10*  to  prevent  unwanted  low 
effects.  (c)  Linear  rotor  speeds  comparable  with  the  axial  gap 
flow  velocity,  i.e.,  of  the  order  of  100  m/sec.  (d)  Provision  for 
controllable  inlet  flow  swirl  with  tangential  velocities  of  the  same 
order  as  the  rotor  speed,  or  the  axial  gap  velocity.  Secondary 
specifications  include  modularity  (i.e.,  a  capability  to  include 
various  types  of  seals  and  seal  lands,  as  well  as  a  range  of  number  of 
glands  per  seal),  a  reasonably  large  nominal  gap  (to  allow  precise 
setting  of  the  eccentricity)  and  solid  construction  to  minimize 
vibrations. 

The  test  section  will  operate  in  air,  discharging  into  the  room. 
For  steady  operation  we  will  directly  connect  to  an  existing  Gas 
Turbine  Laboratory  oil-free  air  compressor  (0.3  kg/sec,  6  atm).  For 
tests  requiring  higher  flow  we  can  operate  in  a  blow-down  mode  using 
the  60  m3  of  pressurized  air  (6  atm)  available.  For  example,  for  4  atm 
inlet  pressure,  a  20  sec  blow-down  will  only  vary  inlet  conditions  by 
about  5 X.  Table  1  lists  several  test  conditions  achievable  with  this 


rig: 
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TABLE  1  Test  Conditions  for  Seals  Test  Rig 


P  /? 
o  e 


<5 

Nominal 


- - -  nAiai 

(No .of  strips)  (Press. ratio)  (atm)  gap  (mm.)  (Kg/sec)  gap  R  No 


Axial  Operating 
P  R.  No.  Mode 


2 

1.54 

1.54 

0.6 

0.127 

15,000 

Steady 

2 

1.54 

1.54 

0.3 

0.169 

20,000 

Steady 

2 

1.54 

1.54 

1.2 

0.254 

30,000 

Blow- Down 

20 

3.85 

3.85 

0.6 

0.127 

15,000 

Steady 

20 

3.85 

3.85 

0.8 

0.169 

20,000 

Steady 

20 

3.85 

3.85 

1.2 

0.254 

30,000 

Blow-Down 

A  conceptual  design  including  a  set  of  nested  eccentric  bearings  for 
whirl  generation  was  included  in  our  Proposal  for  renewal  of  this 
Grant.  That  design  has  since  evolved,  although  retaining  the  same 
operating  principles.  Fig.  1  shows  our  current  preliminary  design. 
The  spinning  shaft  is  supported  by  a  set  of  inner  bearings,  whose  outer 
race  rests  on  an  eccentric  insert.  This  insert  is  housed  on  the  inside 
of  an  intermediate  cylinder,  which  can  itself  spin  concentrically  on  a 
set  of  outer  bearings  supported  by  the  outer  casing.  The  intermediate 
cylinder  rotation  is  provided  by  a  separately  driven  pulley.  This 
arrangement  forces  the  geometrical  axis  of  the  shaft  to  execute  a  small 
satelliting  motion  about  the  casing  center,  at  the  frequency  of  the 
intermediate  cylinder  rotation.  Superimposed  on  this  whirl,  there  is 
the  shaft  spin,  driven  through  a  flexible  connector  that  can  accomodate 
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the  small  parallel  translations  of  the  shaft.  The  linear  speed  of  100 
m/sec  at  the  test  seal  requires,  with  the  30  cm  diameter,  a  spin  of 
6370  RPM.  Since  whirl  occurs  typically  at  1/2  the  shaft  rotational 
speed,  we  plan  to  provide  belt  drive  speeds  up  to  3180  RPM. 

The  air  supply  head  is  structurally  separated  from  the  rotating 
assembly  casing,  to  minimize  vibration  and  stress  transmission.  The 
rotating  assembly  is  supported  on  a  single  ring  flange  for  the  same 
reasons.  A  dummy  disc  is  mounted  on  the  shaft  end  opposite  the  test 
seal  disk  to  provide  mass  symmetry.  To  reduce  vibratory  loads  due  to 
the  eccentricity  of  the  shaft,  we  will  add  counterweights  attached  at 
appropriate  sites  on  the  outer  surface  of  the  intermediate  cylinder. 

Air  is  supplied  to  the  test  seal  inlet  by  two  different  routes, 
one  arranged  to  feed  concentrically  the  seal  inlet  annular  plenum,  the 
other  arranged  to  provide  tangential  momentum  to  that  plenum.  By 
varying  the  relative  feed  rates,  inlet  swirl  can  be  controlled  over  a 
wide  range.  High  pressure  air  is  also  diverted  to  the  back  side  of  the 
test  seal  disk  in  order  to  provide  at  least  partial  thrust  balancing. 

The  seal  land  is  a  separately  mounted  ring  which  can  easily  be 
interchanged  (as  can  the  test  seal  itself).  The  primary 
instrumentation  will  consist  of  two  or  more  high  response,  flush 
mounted,  pressure  transducers  per  seal  chamber .  To  increase 
sensitivity,  these  transducers  will  be  referenced  to  a  ring  manifold 
kept  at  the  average  chamber  pressure  by  means  of  a  number  of  small 
orifices  distributed  over  the  seal’s  periphery  and  connecting  to  the 
manifold.  Additional  instrumentation  will  include  precise  displacement 
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transducers  and  case-mounted  accelerometers,  plus  whirl  and  spin 
tachometers,  and  flow  meters. 

(b)  Analytical  Studies 

The  work  presented  in  the  Appendix  to  our  previous  yearly  report 
has  been  substantially  extended  by  a  careful  examination  of  the 
physical  mechanisms  responsible  for  the  computed  seal  forces.  This  was 
accomplished  by  non-dimens ionalizing  the  governing  linearized  equations 
and  comparing  the  relative  magnitudes  of  the  various  terms.  We  could 
thus  identify  two  primary  fluid  effects  as  being  responsible  for  the 
bulk  of  the  cross-forces:  (a)  Unsymmetric  tangential  momentum 
addition  to  each  chamber  when  the  flow  swirl  varies  axially  in  the 
presence  of  shaft  offset  and  (b)  Azimuthal  variations  of  flow  area 
within  each  chamber  due  to  the  offset.  Both  of  these  induce  pressure 
non-uniformities  that  integrate  to  yield  cross-forces.  The  direct 
contribution  of  shear  forces  is  shown  to  be  small,  but  friction  does 
play  an  important  role  in  generating  the  flow  swirl  changes  that  lead 
to  effect  (a)  above.  The  Alford  passage-contraction  effect  is  seen  to 
contribute  only  to  direct  stiffness,  and  is  thus  of  only  minor 
importance  to  stability. 

The  revised  version  of  our  1984  paper  (Ref.  2),  containing  the 
analysis  that  leads  to  the  above  conclusions,  is  included  as  Appendix  A 
to  this  Report.  This  paper  is  being  submitted  for  publication  in  the 
ASME  Journal  of  Engineering  for  Gas  Turbines  and  Power.  The  initial 
formulation  and  the  numerical  results  (about  50%  of  the  paper)  is 
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identical  with  Ref.  2. 

Also  completed  during  this  work  period  and  listed  as  Ref.  3  is  the 
Master's  Thesis  by  Otto.  W.K.  Lee,  which  contains  the  development  and 
initial  numerical  results  of  the  linearized  theory.  This  was  submitted 
in  Feb.  of  1984.  The  code  was  later  refined  and  the  computations 
repeated  by  a  new  graduate  student,  Eva  Czajskowski,  who  also  performed 
supporting  experimental  studies  on  the  details  of  the  intra-cavity 
pressure  distribution  in  a  seal.  This  work  will  be  reported  in  our 
next  Yearly  Report. 


2.2  Structual  Dynamics 

A  structural  dynamic  model  of  rotating  turbomachinery  has  been 
developed  to  study  the  effects  of  labyrinth  seal  forces  on  dynamic 
stability.  The  model  consisted  of  a  flexible  rotor  with  bearings  and 
labyrinth  seals  mounted  on  a  rigid  or  flexible  case.  In  addition  to 
the  conventional  "stiffness"  and  "damping"  expressions,  additional 
"mass  inertia"  and  "dynamic  memory"  effects  are  used  to  represent  the 
labyrinth  seal  forces.  These  additional  expressions  were  obtained  by 
fitting  the  conventional  frequency  dependent  stiffness  and  damping 
terms  with  Pade'  approximants,  to  give  frequency  independent  (constant 
coefficient)  stiffness,  damping,  mass  and  lag  terms  which  can  then  be 
used  in  a  traditional  eigenvalue  analysis  to  examine  possible 
rotordynamic  instabilities. 

The  above  rotordynamic  formulation  was  developed  in  a  Master’s 
thesis  by  Y.P.  Chang  (Ref.  4),  and  applied  to  the  High  Pressure  Fuel 
Turbopump  (HPFTP)  of  the  Space  Shuttle  main  engine.  Preliminary 
results  indicated  that  the  labyrinth  seal  forces  induce  instability, 
but  this  instability  can  be  delayed  by  increasing  bearing  direct 
damping. 

In  the  analysis  of  Ref.  4,  the  rotor  was  modeled  as  a  rotating 
free-free  beam,  with  rigid  body  translations  q*  ,  q,  and  rotations  B„ , 
8t  as  shown  in  Fig.  2.  The  bending  deflections  W,  V  of  the  rotor  were 
first  described  with  reference  to  a  rotating  frame,  and  later 
transformed  to  a  fixed  frame.  The  rotor  was  interconnected  to  a  rigid 


or  flexible  case  at  the  bearings  and  at  the  labyrinth  seals  by 


appropriate  spring,  damping,  and  seal  forces.  Figure  3  shows  the 
dynamic  model  of  the  High  Pressure  Fuel  Turbopump  (HPFTP)  rotor  used 
for  the  analysis  with  the  associated  bearing  and  seal  locations.  This 
rotor  model  as  well  as  a  model  of  the  surrounding  case,  are  given  by 
Muller  (Ref.  5),  together  with  the  appropriate  dynamic  properties. 

The  labyrinth  seal  forces  for  the  impeller  seals  1,2,3  were 
estimated  using  the  theory  developed  earlier  by  Lee,  Martinez-Sanchez 
and  Czajkowski  (Ref.  2).  In  this  theory,  the  forces  in  the  y  and  z 
directions  are  given  as 

Fy  —  ky  y  y  —  by  y  y  —  k2  b2  y  z 

(1) 

Fa  =k2yy  +  bjyy  —  kyyZ—byyZ 

where  circular  symmetric  motion  of  a  symmetric  shaft  has  been  assumed. 
The  conventional  spring  and  damping  coefficients  kyy,  kzy,  byy,  bzy  are 
dependent  on  the  whirl  frequency  w  as  well  as  the  rotation  frequency  Q. 
Since  the  whirl  frequency  w  is  generally  to  be  determined  from  the 
roots  of  the  resulting  stability  analysis,  an  iteration  process  must  be 
done  to  match  the  roots  with  the  w  assumed  for  the  seal  force 
coefficients.  To  avoid  this  iteration,  the  seal  force  coefficients 
were  fitted  by  Pade  approximants  over  a  wide  range  of  whirl  frequencies 
w,  as  is  often  done  in  aircraft  aeroelastic  calculations  (Ref.  6).  The 
resulting  seal  forces  of  Eq.  (1)  were  then  alternately  expressed  as, 
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where  all  coefficients  are  now  independent  of  whirl  frequency  w,  and 
additional  "mass"  and  "lag"  terms  appear  in  addition  to  the 
conventional  stiffness  and  damping  terms.  Figure  4  shows  a  typical 

example  of  the  2nd  impeller  seal  coefficients  kyy,  byy . and  their 

variation  with  whirl  frequency  w  at  the  rotational  speed  (2  =  3708 
rad/sec,  along  with  the  appropriate  Pade  fit  by  the  whirl  frequency 

independent  coefficients  klyy,  btyy, . 

The  whirl  frequency  independent  seal  forces  of  Eq.  (2)  were 
combined  with  the  structural  dynamic  analysis  to  form  a  conventional 
linear  system  of  constant  coefficient  equations  of  the  form, 

M*Q  +  Bq  +  Kq  +  GW  =  F(t) 

<0  <v  '■O  ^  »v 

Nw+Hw=-Rq  (3) 

<V  <V  /V  /w  4 

where  3  represents  the  structural  degrees  of  freedom,  and  w  are 
augmented  state  variables  defined  by, 
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rt  -fci  (t  -  r) 

w1  (t)  =  e  qi  (r)dr  (4) 


To  examine  Eqs.(3)  for  stability,  the  forcing  terms  F(t)  were  set  to 
zero,  and  the  equations  were  rewritten  in  first  order  form  as, 


where, 


(5) 


0 

I 

0 

3 

-  K 

A# 

-M" 1  B 

-M" 1 G 

<V  M 

,  x  = 

*  M 

• 

3 

0 

N*  1R 

iv  <V 

-N-tH 

IV  »v 

w 

Assuming  homogeneous  solutions  X(t)  =  XePt  led  to  a  standard  eigenvalue 
problem  for  the  roots  p,  which  was  solved  at  each  rotational  speed  ft. 
Any  positive  root  or  root  with  a  positive  real  part  indicated 
instability. 

Figure  5  shows  a  root  locus  plot  of  the  roots  p  for  the  rotor 
mounted  on  a  rigid  case.  Instability  is  seen  to  occur  at  a  rotation 
speed  below  ft  =  1000  rad/sec  and  with  a  whirl  frequency  w  =  2100 
rad/sec.  The  addition  of  some  bearing  damping  is  seen  to  stabilize  the 
system  somewhat.  Figure  6  shows  the  corresponding  roots  for  the  rotor 
on  a  flexible  case.  The  numerical  results  for  these  analyses  seem  to 
show  some  minor  discrpancies  and  discontinuities,  which  will  be 
examined  further  in  a  subsequent  report  by  another  graduate  student. 
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3.  ACCOMPLISHMENTS 

1.  A  labyrinth  seals  test  rig  has  been  designed,  and  is  currently 
being  constructed.  The  test  rig  will  provide  realistic  pressure  drops 
and  Reynolds  numbers,  and  will  attempt  to  measure  both  stiffness  and 
damping  forces  for  various  whirl  and  rotation  frequencies. 

2.  An  analytic  model  has  been  developed  for  labyrinth  seal  forces  in 
turbomachinery.  A  Master’s  Thesis  based  on  this  work  has  been  written 
by  Lee  (Ref.  3)  and  a  paper  was  presented  to  the  3rd  Workshop  on  Rotor 
Dynamic  Instability  Problems  in  High  Performance  Turbomachinery  (Ref. 
2) .  Subsequent  extension  of  this  work  including  non  dimensionalization 
of  the  resulting  equations  has  revealed  the  importance  of  unsymmetrical 
tangential  momentum  addition  to  each  chamber,  and  azimuthal  variations 
of  flow  area  within  each  chamber  due  to  shaft  offset.  A  paper  based  on 
this  work  is  being  submitted  for  publication  to  the  ASME  Journal  of 
Engineering  for  Gas  Turbines  and  Power,  and  is  included  as  Appendix  A 
of  the  present  report. 

3.  A  structural  dynamic  stability  formulation  has  been  developed  for  a 
rotor  system  using  the  labyrinth  seal  forces  developed  earlier  for  the 
seals.  The  whirl  frequency  dependent  seal  forces  were  reformulated 
using  Pade  approximants  as  frequency  independent  stiffness,  damping, 
mass  and  lag  coefficients  which  can  then  be  used  in  a  traditional 
eigenvalue  analysis  to  examine  possible  rotor  dynamic  instabilities. 
The  dynamic  stability  formulation  was  applied  to  the  High  Pressure  Fuel 


1 


I 


Turbopunp  of  the  Space  Shuttle  main  engine.  Further  details  of  this 
work  are  given^in  a  Master's  thesis  by  Chang  (Ref.  4). 
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Fig.  1  Preliminary  Seals  Test  Rig  Design 


Fig.  3  High  Pressure  Fuel  Turbopump  (HPFTP)  Rotor  Model 
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A  linear  analytical  model  for  the  prediction  of  fluid  forces  in 
labyrinth  seals  is  presented  and  discussed.  Numerical  comparisons  to 
literature  data  on  the  important  corss-stiffness  ky,  shows  agreement  to 
within  20 X.  Two  dominant  fluid  effects  are  identified:  (a) 
Unsymmetric  tangential  momentum  addition  to  each  chamber  when  the  flow 
swirl  varies  axially  in  the  presence  of  shaft  offset,  and  (b) 
Azimuthal  variations  of  flow  area  within  each  chamber  due  to  the 
offset.  Both  of  these  induce  pressure  nonuniformities  that  integrate 
to  yield  cross-forces.  The  direct  contribution  of  shear  forces  is 


shown  to  be  small,  but  friction  does  play  an  important  role  in 
generating  the  swirl  changes  that  lead  to  effect  (a)  above.  The  Alford 
passage-contraction  effect  is  shown  to  contribute  only  to  direct 


L.  Introduction 


Fluid-dynamic  forces  arising  from  nonuniform  pressure  patterns  in 
labyrinth  seal  glands  are  known  to  be  potentially  destabilizing  in  high 
power  turbomachinery.  A  well  documented  case  in  point  is  that  of  the 
space  Shuttle  Main  Engine  turbopumps  (ref.  1),  and  other  examples  can 
be  found  in  the  literature,  as  for  instance  in  the  reviews  of  Ehrich 
(ref.  2)  and  Childs  and  Ehrich  (ref.  3).  Seal  forces  are  also  an 
important  factor  for  the  stability  of  shrouded  turbines,  acting  in  that 
case  in  conjunction  with  the  effects  of  blade-tip  clearance  variations 
(ref.  4,5,6). 

The  basic  mechanisms  which  produce  the  uneven  pressure  distrubtion 
in  a  labyrinth  have  been  qualitatively  or  semi-qualitatively  discussed 
in  many  early  references  (Ref.  7,  for  instance).  In  most  instances, 
the  basic  agent  is  found  to  be  flow  swirl  in  the  glands,  either  from 
?re-swirl  (as  in  the  case  of  turbine  shrouds)  or  from  frictional 
interaction  with  the  rotating  shaft  (as  in  multichamber  jet  engine 
seals).  Quantitative  modeling  of  these  forces  has  also  been  reported 
by  a  number  of  authors  (refs.  8,  9,  10),  using  lumped-parameter  models 
for  each  gland.  These  models  yield  in  general  predictions  of  the 
direct  and  cross-wise  stiffnesses  and  damping  coefficients  for  small 
shaft  displacements,  and  are  useful  for  linear  stability  analyses. 
Non-linear  predictions  for  fully  developed  unstable  operation  are  less 
advanced . 

Fairly  extensive  data  also  exist  on  the  stiffness  factors  of  seals 
of  various  geometries  (refs.  4,  5,  11,  12,  13).  These  have  been 
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generally  obtained  in  rotary  rigs  with  adjustable  shaft  eccentricity. 
Much  less  satis-factory  is  the  situation  with  respect  to  data  on  damping 
coefficients  due  to  labyrinth  seals,  since  these  require  dynamic 
measurements  on  either  vibrating  shafts,  or  shafts  fitted  with 
adjustable  whirl  mechanisms.  Yet  these  data  are  almost  as  essential  as 
those  on  stiffnesses,  since  the  corresponding  induced  forces  are  of  the 
same  order.  Ref.  (13)  reports  damping  data  for  a  labyrinth,  but  with 
no  positive  control  or  measurement  of  inlet  whirl. 

In  this  paper  we  report  on  the  development  of  a  linear  model  for 
the  prediction  of  labyrinth  seal  forces  and  on  its  comparison  to 
available  stiffness  data.  We  also  present  a  discussion  of  the 
relevance  of  various  physical  effects  that  lead  to  fluid  forces  on  the 
labyrinth. 

2.  Model  Formulation 

The  model  is  very  similar  in  its  main  outline  to  those  of  Kostyuk 
(ref.  8)  and  Iwatsubo  (ref.  9).  It  describes  the  flow  of  an  ideal  gas 
through  the  seal  chambers,  assuming  largely  constant  temperature,  but 
allowing  for  isentropic  acceleration  towards  the  narrow  gaps  and  also 
for  isentropic  azimuthal  flow  redistribution  in  each  chamber.  Each 
chamber  is  assigned  a  pressure  P;  and  azimuthal  velocity  C! ,  and  these 
quantities  are  governed  by  equations  of  mass  and  azimuthal  momentum 
conservation,  written  in  integral  form. 

The  axial  flow  rate  q  through  each  seal  throttling  is  approximated 
by  a  commonly  used  expression,  which  basically  derives  from  Bernoulli’s 


equation  with  a  density  halfway  between  those  of  the  two  adjacent 
chambers.  Per-unit  length  in  the  azimuthal  direction,  this  gives 


qj  =  ^ 


where  5,  is  the  width  of  the  narrow  gap  (see  fig.  1)  and  =  c,.  ft  is 

i 

the  product  of  the  usual  contraction  coefficient  c*.  times  a  carryover 

i 

factor  pt  to  account  for  nonzero  upstream  axial  velocity  and  nonzero 
pressure  recovery  in  the  downstream  chamber.  Eq.  (1)  is  assumed  to 
apply  locally  at  each  time  t  and  azimuth  This  semi-incompressible 
approximation  is  known  to  be  reasonalbe  up  to  gap  Mach  numbers  of  about 
0.5;  however,  the  last  gap  or  two  of  a  labyrinth  with  a  high  overall 
pressure  ratio  may  be  above  that  Mach  number,  and,  in  particular,  the 
last  chamber  may  choke.  We  have  partially  accommodated  this  effect  by 
retaining  Eq.  (1)  throughout,  but  replacing  it  by  a  choked-flow 
expression  in  the  last  chamber  only  if  the  first  approximation 
indicated  sonic  or  supersonic  conditions  there. 

With  reference  to  the  geometry  of  fig.  1,  the  governing  equations 
within  each  chamber  are  those  of  mass  and  azimuthal  momentum 
conservation: 


ft (pi  f* }  +  iw (P|  f|  c* )  +  q  ♦ 1  ~  q‘  =  0 


|r  (Pi  f(  c,  )  +  f-  (Pi  f;  c,  )  +  q;  ,  i  C,  -  q,  C|  _ !  +  r,  "IT 


These  equations  are  first  linearized  about  a  condition  of  zero 
eccentricity.  The  zero’th  order  approximation  provides  a  basic  flow 

f  *  * 

rate  q  and  pressure  and  azimuthal  velocity  distributions  Pj  ,  c,  , 
(Appendix  1).  The  first  approximation  then  provides  linear  equations 
for  the  perturbations,  defined  by 

Pi  =  pf  (l  +  Ki )  ;  qi  =  qf  (l  +  $i )  ;  c,  =  c*  (l  +  m  )  (4) 

where  ,  m  are  functions  of  t  and  w  =  Rs  The  right  hand  sides 

of  these  equations  are  determined  by  an  assumed  accentric  motion  of  the 
shaft,  whose  center  follows  a  circular  path 

,  .  iflt 

x,.  +■  iyc  =  r  e  (5) 

where  Q  is  the  shaft  vibration  frequency,  closely  identified  with  one 
of  its  natural  frequencies.  This  results  in  a  gap  distribution  (see 
Fig.  2) 

5i  =  6*  -  r  cos(#  -  Qt)  (6) 

where  ^  -  Qt  is  the  angle  measured  in  the  whirl  direction  from  the 
narrowest  gap.  The  perturbation  quantities  ,  m  and  then  have 
similar  sinusoidal  variations  in  time  and  space,  with  whirl-dependent 
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amplitudes  and  phases: 

=  Re[t,e_l(*“0t)],  etc  (7) 

/\ 

where  Re()  indicates  real  part  and  is  a  complex  amplitude  factor. 

The  process  of  linearizing  Eqs.  (1),(2)  and  (3)  is 
straightforward,  but  tedious  (see  Refs.  14,15).  After  eliminating  flow 
variations  ($|  )  in  favor  of  pressure  ($j  )  and  tangential  velocity  (rji  ) 
variations,  one  obtains  a  system  of  two  linear  equations  per  labyrinth 
gland  as  follows: 


riF, 

0  1  l 

-Bj  +  QA]  **  iEj 

~Ci  1 

liQi 

is‘J  .r 

.  -Mi  +  OK*  -  iP* 

-N,  +  flL,  -  iR,  J 

Hi 


f  iDi  01  Mim  If  -Z4  iG,  ) 

l  iO,  0j  \  *  I  '  \  ~Y,  +  nx,  -  Wt  f 


(8) 


Hi  <• ; 


The  expressions  for  the  coefficients  A;  ,  Bj  . ,  are  given  in 

a  A 

Appendix  2.  For  the  first  chamber  (I  =  1),  - 1  ,  rji-i  represent  the 
inlet  non-uniformity,  and  for  the  last  chamber  (i  =  k),  they  represent 
the  outlet  nonuniformity.  For  the  calculations  reported  here,  these 
were  taken  to  be  zero. 

If  the  pressure  perturbations  Pf $|  (4,t)  are  calculated  for  each  of 
the  gaps,  then  the  net  pressure  forces  are  obtainable  by  integration 
around  the  rotor  periphery,  followed  by  summation  for  all  chambers.  It 
is  advantageous  to  project  the  forces  in  the  directions  towards  the 
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o 

instantaneous  ainisua  gap  (Fd  )  and  90  ahead  of  it  in  the  whirl 
direction  (F, )T  as  shown  in  Fig.  2.  We  obtain  for  each  chamber 

Fd  +  F,  =  -R,^  Pr?,el(#"nt)d^  =  -nRt«lpf?i  (9) 

i  i  ° 

The  forces  F, ,  F„  along  the  fixed  coordinate  directions  are  then 
given  simply  by 

F,  +  iF,  =  (Fd  +-  iF,  )  elfit  (10) 

till 


3.  The  Stiffness  and  Damping  Coefficients 

The  most  general  presentation  of  the  results  of  a  linear  theory 

such  as  ours  would  be  plots  of  the  direct  and  quadrature  forces  per 

F  F 

unit  whirl  amplitude  versus  whirl  speed  12,  with  rotation  speed  w 

and  other  variables  as  parameters.  Polynomial  approximations  for  Fd (0) 
and  F, (fl)  are  useful  in  condensing  this  information.  The  results 
presented  here  have  been  cast  in  the  form  of  a  linear  approximation  in 
(2,  leading  to  the  definition  of  stiffness  and  damping  coefficients 
according  to  the  scheme 

I  F,  =  -  Kf*Xc  +  Kiyyc  -  Cfl  Xc  +  Cx,yc  (11) 

'  F»  -  ~  Kx  y  Xc  _  Kj  i  yc  ~  ,  Xc  ~  C, ,  yc 

or,  equivalently,  using  Eqs.  (5)  and  (10), 
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=  _(KlI  +  iKf  *  )  -  ifl(Clf  +  iCiy  )  (12) 

Since  such  a  representation  is  a  truncation  of  some  more  complex 
behavior,  the  K  and  C  coefficients  must  be  regarded  as  functions  of  (2. 
However,  our  results  indicated  fairly  small  deviations  from  linearity 
in  (2  (see  Appendix  3).  Also,  the  linear  approximation  may  be  local 
(tangent)  or  global  (secant),  as  illustrated  schematically  in  Fig.  3. 
In  the  results  presented  here,  the  "secant"  approximation  has  been 
used: 


K  =  .  Ldm-Ll±k£l  c  -  f,(_Q)  -  f,  (~n) 

RII  Or-  90r 


***  2r  C,»  2flr 


(13) 


This  can  be  shown  to  be  equivalent  to  a  commonly  used  procedure 
(Refs.  (8),  (9),  (15)  in  which  the  equations  are  solved  for  a  linearly 
vibrating  shaft,  and  the  calculated  forces  are  separated  into  those  in 
time  phase  with  the  displacement  (K’s)  and  those  at  90°  (C’s). 

A  very  simplified  dynamic  model  will  help  to  put  in  perspective 
the  roles  of  the  different  coefficients  in  stability  analysis. 
Assuming  a  shaft  with  mass  M  and  structural  stiffness  K<,  ,  the  equations 
of  motion  for  small  side  displacements  x,  y  can  be  combined  into 


Mz  +  Cz  +  Kz  =  0 
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A  A 

where  C  =  CfI  +  i  Cx  y  ,  K  =  Ko  +  Kx  x  +  i  Kx  y  and  z  =  x  +  iy. 


Assuming  K,  x  ,  Kjj,,  C|X/4M  and  C\  y  /4M  are  all  small  compared  to  Kj  ,  as 
is  likely  to  be  the  case  in  practice,  we  can  define  the  (small) 
nondimens ional  parameters 


=  ILl 
Ko 


Ski 


= 

2v/kTm 


Then  a  simple  analysis  shows  that,  to  the  first  approximation,  the 

0  t 

shaft  complex  displacement  z  will  vary  as  e  9  ,  where 


-Ztl  -  (  -  Cix  *  I  *  i  1  *  i  kxf  -  $ry) 
v^T/M  *  * 

Thus  any  nonzero  kxy  will  be  destabilizing  (in  one  or  the  other 
whirl  direction),  while  a  negative  $xx  will  be  always  destabilizing; 
kxx  and  $xy  will  simply  modify  the  shaft  natural  frequency.  Also,  the 
effects  of  equal  values  of  |Kxy  |  and  -(2  Cxx  are  seen  to  be  equivalent. 

This  discussion  has  served  to  indicate  that  knowledge  of  the 
damping  factor  Cxx  is  at  least  as  essential  to  studies  of  fluid-induced 
destabilizing  forces  as  is  the  side  force  factor  K,  y .  Yet,  due  to  the 
more  difficult  experimental  conditions,  much  fewer  data  are  available 
on  Cx  x  than  on  Kx  y . 


L.-- 
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4.  Discussion  of  the  Physical  Effects  Involved 


4.1  Non-Dimensional  Formulation 

For  purposes  of  analysis  it  is  advantageous  to  recast  Eqs.  (8)  in 

non-dimensional  form.  It  is  also  convenient  to  replace  the  second  of 

the  two  equations  in  (8)  (i.e. ,  the  momentum  equation)  by  the 

combination  resulting  from  subtracting  from  it  the  first  equation 

* 

(continuity),  times  q  . 

the  list  of  dimensionless  parameters  used  contains  a  few  purely 
geometrical  ratios,  namely 


(14) 


plus  a  number  of  flow-related  parameters.  These  are: 


(a)  A  parameter  characterizing  the  gland-to-gland  pressure  ratio,  or 
the  restriction  axial  Mach  number: 


(15) 


If  M,  is  the  axial  Mach  number  'defined  with  the  chamber  speed  of 


;v  V'v - 


sound,  -ARf  T  ),  in  the  (i-1,  i)  restriction,  one  can  show  that,  within 
the  model  appreximations, 


q 


Jr  M, 


(16) 


(b)  A  parameter  related  to  the  flow  swirl  in  each  gland,  or  to  the 
angle  between  the  flow  in  the  restriction  and  the  axial  direction: 


a.  = 


*  * 
BlLSl. 

q* 


1  * 

T/ij  MI 


(17) 


where  Mt  is  the  tangential  Mach  number,  Mt  =  C;  /  /rR*T 


(c)  A  mixed  parameter  which  compares  pressure  gradients  to  dynamic 
head  gradients: 


Hi  = 


*  * 

q  Ci 


i-z* 


M,  Mt 


(18) 


(d)  A  "flow  twisting"  parameter  which  measures  the  axial  swirl 


gradients: 


Finally,  the  spin  and  whirl  angular  velocities  of  the  shaft  are 


expressed  by  the  parameters 


C| 


(20) 


In  terms  of  these  variables,  the  set  of  equations  (8)  becomes: 


,,  .  1  4  t  .  r  1+5?  .  .  Li  ,  :  0?  -  1  ?  .  L, 

-d  +  “ )  Ci  -  x  +  [-ZT"  +  i  a{  (w,  -  1)  ]  Ci  -  ^ -  Cm  i  -  i  -J-  o,  q, 

q*  ql  ‘  ql  • 


=  s,  [a,  +  i  Li  a{  (w,  -  1)] 


(21) 


*  <l  *  <  -  *■  *  stht  [l  - ~ 1)!  i  - iff'  £-><■ 


a  -  g, )  i,.,  (-  -  i)] 


*■  i  T1-  a,  (wt  -  1)>  fft  =  * i  gi 


(22) 


i 

In  order  to  assess  the  numerical  range  of  the  various  terms  in 
these  equations,  it  is  useful  to  particularize  to  the  case  of  a  simple 
1-cavity  (two  restrictions)  seal.  If  the  upstream  and  downstream  fluid 

/\  /\  'N 

variables  are  assumed  uniform,  then  Ci-i  =  Cmi  =  Hi  - 1  =  0  in  this 
case.  In  what  follows  we  will  drop  the  subscript  i,  since  only  the 
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value  i=l  is  involved.  This  case  also  allows  a  relatively  simple 

*  * 

calculation  of  the  chamber  swirl  c  (and  hence  of  g  =  1  -  ^  ):  the 

c 

steady-state  momentum  equation  (Eq.  A-3)  gives  a  zerolh  order  solution 
$  I 

c  —  cQ  when  friction  is  neglected.  Iterating  with  this  value  then 
gives  an  improved  estimate  as  follows: 


g  =  1 


* 

_  So.  - 
* 
c 


A'2a. 
8  M 


[1  ~ 


A"  U" 
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(So  -  1)*] 
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4.2  Orders  of  Magnitude  and  Simplified  Equations 

For  concreteness,  we  will  now  refer  to  the  results  reported  in 
Appendix  3  for  a  particular  choice  of  parameters.  We  find  in  this  case 
friction  factors  A'  ~  A'*  -  0.006,  and  axial  and  tangential  Mach 

numbers  ^  -  0.61,  Mt  ^  0.28  (at  an  inlet  swirl  velocity  c^  =  100 
m/sec).  This,  plus  the  dimensions  chosen,  gives: 

A  =  0.083  L  =  0.053  H  =  0.38  a  =  0  P  =  l 
q  =  0.83  a  =  0.28  U  =  4.4  g  =  -0.025  w  =  0.75 


Using  these  values  it  is  found  that  the  terms  involving  the 
friction  factors  A',  A '  •  in  Eq.  (22)  are  about  an  order  of  magnitude 
smaller  than  other  competing  terms,  and  can  be  neglected,  with  the 
important  exception  that  the  whole  RHS  of  Eq.  (22)  is  proportional  to 
g,  and  hence  to  the  friction  factors  (Eq.  23),  and  cannot  a-priori  be 
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ignored.  Eqs.  (21)  and  (22)  reduce  to: 


f1  +J*'  +  i  j  °(w  ~  D] 


A  Ii  ^ 

<  -  i  —  a  rj  -  e{a  +  iLo(w  -  1)] 


il7  2?  +  [l  +  i|  o(w  “  1 )  ]  h  -  eg 


(24) 


(25) 


This  can  yet  be  reduced  by  noting  that  70  (w  -  1)  ^  -  0.05;  if  this 

A 

tens  is  dropped  from  the  bracketed  terms,  we  obtain: 


£  _  €  a  +  i  Lo(w  -  1  +  a/A) 
1  +  &  .  _  L* 

— - e  +  /70  TT 

q2 


A1 


n 


e 


1  ♦  &■ 

_JL_ 


g  -  Ha  7— (w  -  1)  +  iaJT  7 
_ A _  A 


1J_^  + 

q2 


(26) 


(27) 


4.3  Interpretation  of  Physical  Effects 

A 

Since  F,  +  iF,  is  proportional  to  -(  (Eq.  9),  we  see  from  Eq.  (26) 
that  the  "passage  contraction  factor",  a,  (Eq.  13)  can  only  contribute 
to  Fd ,  which  is  dynamically  much  less  significant  than  the  cross-force 
F,  .  This  result  is  to  be  contrasted  with  the  simplified  model 
postulated  originally  by  Alford  (Ref.  7),  where  axial  decrease  of  flow 
gaps  was  identified  as  a  major  source  of  turbomachine  whirl.  The 
explanation  is  to  be  found  in  the  rapid  pressure  redistribution  by 
sound  waves  within  a  chamber,  which  makes  M,  and  qj  relatively  small 
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and  prevents  the  second  term  in  the  coefficient  of  $  in  Eq.  (24)  from 

A 

contributing  a,  significant  imaginary  part  to  t,  due  to  a  nonzero  a. 
This  conclusion  would,  of  course,  be  reversed  for  cases  where  —  =  — ^ 
is  much  less  than  taken  in  our  example.  Notice  also  that,  within  the 
approximations  adopted,  this  "Alford  seal  effect"  is  the  only 
contributor  to  direct  (displacement-aligned)  forces,  and  hence  to  seal 
stiffness.  Other  contributions  would  arise  from  terms  neglected  in 
arriving  at  Eqs.  (24)  and  (25),  but  it  is  clear  from  this  discussion 
that  such  effects  are  numerically  smaller  than  those  leading  to 
cross-forces. 

Concentrating  now  on  F,  ,  we  can  use  Eqs.  (17)  and  (18)  for  o  and  H 
and  substitute  the  imaginary  part  of  Eq.  (26)  into  Eq.  (9),  to  obtain 


F,  =  nR,  «P*t  J 


1  -  w  - 


l  -|1  <1  ♦<*> 


The  factor  1  -  w  -  g/d  in  Eq.  (28)  identifies  the  two  main  driving 

mechanisms  for  cross-force  generation,  the  part  1  -  w  can  be  traced  to 

the  rate  of  change  of  tangential  flow  area,  as  seen  by  the  swirling 

£ 

flow,  and  it  becomes  zero  when  w  =  1,  i.e.,  when  OR,  =  c  .  This  (2 

dependence  constitutes  the  damping  component  (C,,)  of  the  cross-force. 
The  other  part,  g/d,  is  proportional  to  the  change  in  swirl  from  inlet 
to  chamber;  this  is  an  "ejector  pump"  effect,  and  will  be  discussed 
further  in  what  follows.  Numerically  speaking,  in  our  example,  we  have 
1  -  w  =  0.25  and  g/d  =  -0.30,  so  both  these  effects  contribute  about 


equally.  In  terms  of  direction,  we  notice  that  the  area-change  term 

1  -  w  contributes  a  whirling  force  in  the  direction  of  the  flow  swirl 
♦  t 

c  whenever  ft  <  c  /Rs  ,  while  the  term  -g/d  gives  a  force  in  the  whirl 

direction  (F,  >  0)  whenever  the  flow  swirl  is  highest  (in  the  same 

sense  as  the  whirl)  at  entrance,  and  tends  to  be  reduced  by  friction. 

Looking  at  the  denominator  of  Eq.  (28),  we  see  that  the  term 

1  +  01  is  about  five  times  the  other  part  (L*  /4l  p2  -r) .  It  can  be  shown 

that  neglecting  this  latter  term  is  tantamount  to  neglecting  -in  -  (  in 

A 

Eq.  (25).  Since  this  term  arises  from  the  azimuthal  pressure  gradient 
in  the  momentum  equation,  this  shows  that  the  unsymmetric  injection  of 
excess  momentum  (represented  by  the  right  hand  side  of  Eq.  (25)  results 

A 

mostly  in  an  uneven  azimuthal  velocity  (^),  with  a  smaller  effect  on 

/A 

azimuthal  pressure  redistribution  ($)• 

Interpretation  of  the  dominant  effects  is  facilitated  by  restoring 
the  simplified  Eqs.  (24)  and  (25)  to  the  original,  physical  variables. 
Using  primed  symbols  to  represent  the  fluctuating  parts  of  quantities, 
we  obtain 


*  *,df' 

p  c  (5T 


,  *  0t\ 

+  c  - — ) 

dw 


+  pV  —  +  q-  -  q1 

dw  exit 
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(30) 


where  is  that  part  of  the  inflow  variation  which  is  proportional  to 


gap  variation  (the  other  part  is  due  to  the  variations  of  cavity 


pressure,  and  plays  a  smaller  role  in  the  momentum  equation).  If  in 
this  momentum  -equation  (Eq.  (30))  we  ignore  f*  we  obtain  simply 

<7W 

the  statement  that  the  quantity  q(c  -  q,  )  remains  azimuthally  uniform 
despite  the  shaft  eccentricity.  This  means  that  in  places  where  more 
fluid  is  entering  the  cavity  (q1  >  0),  its  velocity  increase  c  -  q, 
gets  reduced,  so  that  no  significant  inertia  force  variation  arises;  to 
the  extent  that  this  balance  is  not  achieved,  (particularly  when  d/L  or 
^  are  small)  azimuthal  pressure  gradients  also  arise.  The  seal  is  in 
some  sense  acting  as  an  azimuthal  ejector  pump,  where  variations  of 
tangential  momentum  addition  result  in  a  combination  of  flow  speed  and 
pressure  gradients. 

With  c'  mainly  determined  from  this  momentum  balance,  Eq.  (29) 
then  indicates  that  the  gradients  of  these  velocity  fluctuations,  plus 
those  of  the  flow  area  (f1),  imply  local  net  mass  inflow  or  outflow 
(q'oiit  ~  9*)-  As  noted  before,  these  have  two  parts,  one  proportional 
to  S',,it  -  S',  which  leads  to  the  Alford  seal  effect,  and  the  other 
due  to  pressure  variations:  if  p'  is  locally  positive,  this  acts  both 
to  inhibit  inflow  and  increase  outflow,  giving  a  net  negative  inflow 
fluctuation.  It  is  this  pressure  fluctuation  that  eventually  produces 
the  bulk  of  the  cross-force  F,  .  Note  also  the  origin  of  the  damping 
component:  the  material  derivative  of  flow  area  (tt—  +  c  - — ) 

<7t  <7W 

$ 

vanishies  if  c  =  flRj ,  and  generally  decreases  linearly  as  n  increases. 

4.4  Shaft  Forces  Due  to  Shear  Stresses 

In  addition  of  the  forces  Fd  and  F,  discussed  so  far,  the  shear 
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stress  fluctuations  also  produce  net  forces  on  the  shaft.  For  the  ith 
chamber, 


(Fd  +  iF,  )t  =  i  U"R,  /a2n  r"*lSd* 


(31) 


Using  our  analysis,  these  can  be  expressed  as 


(Fd  +  iF,), 


A 

i  Ai",~8U~i"  P*C*Z  nR*  (Si_  +  («!“  i) 


(32) 


If  the  Alford  seal  effect  is  ignored  (a  =  0),  our  simplified 

A 

solution  (Eqs.  (26)  and  (27))  show  that,  for  a  one-cavity  seal,  $  is 

A 

imaginary  and  q  is  real  (i.e.,  velocity  fluctuations  are  specially 

aligned  with  the  shaft  displacement,  while  pressure  fluctuations  are 

o  A 

90  out  of  alignment).  This  means  that  only  the  q  (velocity 

fluctuation)  portion  of  Eq.  (32)  contributes  to  (F, )^.  The  ratio  of 

this  to  the  pressure- induced  F,  is 


( F,  )  , ,,  _ 

7=—  =  +  3-  (1  +  2H)M(  (s  -  1)  t  (£*) 

(F,)p  4  K 


(33) 


and,  in  the  numerical  example  of  App.  3,  amounts  to  approximately 
0.0015.  Thus,  frictional  effects  are  secondary  in  directly  determining 
cross-forces,  but  one  should  keep  in  mind  their  role  in  creating  flow 
"twist",  i.e.,  chamber-to-chamber  swirl  variations,  which  do  play  a 
significant  role  in  the  production  of  pressure  forces. 
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5.  Numerical  Results  and  comparison  to  Literature  Data 

5 . 1  Data  Used 

Benckert  and  Wachter  have  published  an  extensive  set  of  data 
(Refs.  11,16)  for  multicharaber  labyrinth  seals  of  simple  "straight"  or 
"full"  types  (Fig.  4).  The  data  were  taken  in  a  static-offset  rig 
operating  on  pressureized  air,  and  induced  forces  were  obtained  by 

integration  of  measured  azimuthal  pressure  variations  on  a  number  of 
seal  cavities.  Labyrinths  with  up  to  23  chambers  were  used.  The 

experiments  allowed  variation  of  shaft  speed,  u,  overall  pressure  ratio 
Patn/Poi  rotor  eccentricity  r,  seal  geometry  (5;,  ,  h;  ,  Fig.  1), 

number  of  chambers  and  entry  swirl  c0  .  The  seal  flow  rate  q  was 

measured  and  an  averaged  carryover  factor  u  was  deduced  from  these  data 
and  reported  in  a  number  of  instances.  In  our  calculations  we  used 
these  "measured"  factors  when  available  directly;  in  other  cases,  we 
adopted  values  measured  for  chambers  of  the  same  geometry,  or,  for  the 
"full"  type  of  seal,  where  little  carryover  is  expected  we  used  ju  =  1. 
In  a  number  of  separate  runs,  we  attempted  an  a-priori  determination  of 
U.  The  contraction  coefficient  ce  was  taken  as  a  function  of  Reynolds 
number  and  strip  geometry  as  given  by  Vermes  (Ref.  17),  and  the 
carry-over  factor  0  was  taken  from  Komotori’s  theory  (ref.  21).  The 
results  of  this  procedure  were  unsatisfactory,  with  the  prediction 
errors  being  several  times  higher  than  when  the  measured  flow 
coefficient  was  used. 


Data  of  Brown  and  Leong  (Ref.  18)  were  also  used  for  validation  of 
the  undisturbed  flow  predictions. 


5.2  Undisturbed  Flow  Parameters 


Figure  5  compares  Brown  and  Leong  s  data  on  the  axial  pressure 
distribution  in  an  11-chamber  test  seal  with  our  calculated  undisturbed 
pressure  distribution.  There  is  good  agreement  except  for  the  sharp 
pressure  drop  shown  by  the  data  between  the  inlet  and  the  first 
chamber.  This  is  probably  a  reflection  of  a  reduced  carry-over  factor 
on  the  first  strip;  the  calculation  used  a  constant  u  (the  value  used 
is  irrelevant  to  the  comparison) . 

Benckert  and  Wachter  reported  for  one  particular  case  the  axial 

variation  of  azimuthal  velocities  c*.  This  was  for  a  23-chamber  seal 
£ 

with  c,  =  40  m/sec,  P8/P„  =  *66,  R,  =  0.15m,  u  =  1000  rad/sec, 
ri  =  0.25  mm,  r*  =0,  5*  =  0.5  mm,  £,  =  4  mm  and  h;  =  6  mm.  The  data 
are  shown  in  figure  6,  together  with  the  code  predictions.  The  good 
agreement  shown  is  important  for  the  prediction  of  disturbance  side 
forces,  which  depend  critically  on  swirl  velocities.  These  results 
appear  to  validate  the  formulation  used  for  the  friction  factors 
between  the  fluid  and  the  stator  and  rotor  surfaces  (turbulent  pipe 
flow  formulae  with  a  standard  correction  for  "pipe"  curvature) 


5.3  Stiffness  Coefficients  Without  Shaft  Rotation 

The  cross-spring  coefficients  =  —  K, y  for  a  number  of  cases 

from  Benckert  and  Wachter’ s  tests  with  a  non-rotating  shaft  were 


i 

! 

calculated  and  the  results  are  summarized  in  Table  2  and  figure  7.  The 
key  in  Table  2>describing  the  test  parameters  is  explained  in  Table  1. 

The  eccentricity  ri  used  in  the  tests  was  0.15  mm,  except  for  Run 
17,  which  had  r,  =  0.25  mm. 

Table  1  Key  for  Table  2  (2nd  column) 

(Type,  h,  6,  «) 


Type 

h :  chamber 
height 

S :  clearance 

<:  pitch 

S  =  Straight-through 

0  =  2.75mm 

0  =  0.25mm 

0  =  5mm 

F  =  Full-interlocking 

1  =  6.25mm 

2  =  0.5mm 

1  — 

3  =  6mm 

2  =  4mm 

As  Figure  7  illustrates,  the  calculated  values  are  somewhat  lower 
than  the  data  (about  19X  for  series  (s, 0,0,0),  5 %  only  for  series 
(s, 1,0,0)).  The  trends  of  the  calculation  are  in  agreement  with  those 
observed  in  the  tests.  In  particular,  Ky,  is  seen  in  figure  8  to  be 
approximately  proportional  to  inlet  swirl  and  to  overall  pressure 
ratio,  both  in  the  tests  and  in  the  calculations  (although,  as 
mentioned,  with  a  somewhat  lower  proportionality  factor  in  the  latter 
case) . 

5.4  Stiffness  Coefficients  with  Shaft  Rotation 

Results  similar  to  those  in  the  previous  section,  but  including 


various  shaft  speeds  are  given  in  Table  3  and  figure  9. 


The 
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eccentricity  is  still  static  (no  whirl,  (2  =  0),  and  is  0.25  mm  in  all 
cases.  The  parameter  Eq  was  used  in  Ref.  11  to  correlate  entry  swirl, 
and  is 
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The  comparison  of  data  and  theory  shown  in  figure  9  indicates  more 
scatter,  but  less  systematic  deviation  than  in  the  cases  without  shaft 
rotation  (figure  7  and  8).  The  agreement  is  best  for  all  the  cases 
with  17  chambers  (solid  symbols  in  figure  9),  which  show  an  average 
error  of  8.5%  and  little  scatter. 

5.5  Discussion 

The  two  principal  sources  of  uncertainty  in  calculations  of  the 
type  reported  here  are  the  friction  factors  (A*,  A")  and  the  carry-over 
coefficient  /J.  The  friction  factor  could  in  principle  be  substantially 
increased  by  the  relative  rotation  of  shaft  and  casing,  since  the  fluid 
in  each  chamber  is  strongly  sheared  and  develops  marked  secondary  flow 


patterns,  leading  to  enhanced  mixing.  Examination  of  data19  for  the 
somewhat  related  case  of  turbulent  pipe  flow  with  swirl  does  indeed 
show  friction  increases  of  up  to  a  factor  of  four  at  high  swirl.  An 
accurate  prediction  of  wall  friction  under  the  complex  flow  conditions 
of  a  labyrinth  gland  is  not  possible  at  this  time,  and  this  is  an  area 
requiring  more  experimental  and  analytical  work. 


The  inpact  of  friction  factor  inaccuracies  on  calculated 
cross-spring  coefficients  could  be  important,  although  not  easily 
general izable.  In  general,  the  cross-forces  increase  with  chamber 
swirl,  and  also,  in  part,  with  swirl  changes  from  chamber  to  chamber 
(terms  Mt  and  g /A  in  Eq.  (28)).  In  a  multi-chamber  seal  with  excess 
inlet  swirl,  the  first  effect  will  make  K*x  decrease  as  the  friction 
factor  is  increased,  since  the  additional  friction  will  accelerate  the 
axial  decay  of  flow  swirl.  On  the  other  hand,  the  changes  in  swirl 
from  chamber  to  chamber  will  be  increased  by  an  increase  in  friction, 
and  this  effect  is  dominant  at  low  friction,  and  especially  for  short 
seals,  where  c  itself  does  not  decrease  much.  An  example  of  this 
behavior  for  a  one-chamber  seal  is  presented  in  Appendix  3,  figure  A2, 
where  it  can  be  seen  that  for  this  particular  case,  increasing  friction 
would  lead  to  increases  in  |K»,|.  The  behavior  typical  of  long  seals 
and  large  friction  is  illustrated  in  figure  10,  corresponding  to  case  1 
of  Table  2;  here,  an  increase  of  the  friction  factor  above  the  nominal 
value  leads  to  a  reduced  Kyx,  although  a  reduction  by  more  than  about 
0.6  would  lead  to  the  same  effect. 

The  carry-over  coefficient  /3  is  clearly  another  source  of 
uncertainty  in  the  model.  Unfortunately,  the  state  of  the  art  in  a 
priori  predictions  of  (3  is  not  satisfactory.  A  fuller  discussion  of 
this  point  is  given  in  reference  24,  which  discusses  the  relationships 
of  the  older  models  of  Vermes17,  Egli.20,  Komotori11  to  more  recent 
models  and  data  (Benvenuti,  et  al.25,  Hendricks2®).  Numerical  methods 
for  internal,  turbulent  and  fully  separated  flows  are  not  yet  capable 


of  better  than  10-20*  accuracy  (see  Ref.  22),  but  full  Navier-Stokes 
computations  are  now  approaching  the  capability  for  modeling  Reynolds 
numbers  of  the  order  of  104  ,  as  required  for  labyrinths.  This  is 
probably  the  best  hope  for  determination  of  A  and  P,  but  for  dynamics 
studies  in  seals,  these  3-D  methods  may,  in  any  case  be  too  laborious, 
and  will  need  to  be  integrated  with  linearized,  overall  models  of  the 
type  described  here. 

6.  Conclusions 

A  linear  analytical  model  for  the  prediction  of  fluid  forces  in 
labyrinth  seals  has  been  presented  and  discussed.  Comparison  to 
literature  data  shows  agreement  to  within  20*  for  the  important 
cross-stiffness  K,,  .  The  dominant  fluid  effects  are  shown  to  be  (a) 
Unsymmetric  tangetial  momentum  addition  to  each  chamber  when  flow  swirl 
varies  axially  in  the  presence  of  shaft  offset,  and  (b)  Azimuthal  flow 
area  variation  due  to  the  offset.  Both  of  these  induce  an  unsymmetric 
pressure  nonuniformity  that  integrates  to  yield  cross-forces.  The 
direct  contribution  of  shear  forces  is  seen  to  be  small,  but  friction 
does  play  an  important  role  in  generating  swirl  changes  that  lead  to 
effect  (a)  above.  Alford’s  passage  contraction  effect  is  shown  to 
contribute  only  to  direct  stiffness. 
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Appendix  1.  The  Unperturbed  Solution 

* 

Squaring  Eq.  (1)  and  adding  for  all  chambers  yields  for  the  nominal 
flow  rate 


—  rj 


Ri T  y  (—hr) 

i.l  Ml  «,** 


(A-l) 


Also,  adding  for  the  first  n  chambers  only  gives 


I  <-ihr> 

i-jjl _ Mi *  4l* 

1  <-i br) 

i  •  i  W 


(A-2) 


The  momentum  equation  (Eq.  (3))  becomes  in  the  steady  state 


with 


q*(cf  -  cf.  t)  +  r,'  U,1 


.  1  *  , ,  *. 
ti  =  -a-  pi  \ v 


r"  =  -i-  pf A."  («R,  -  c*)’ 


r"  U"  = 
Ti  ui 


(A-3) 
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and  A  (the  Darcy  friction  factor)  given  by  a  modified  pipe-flow 
expression 


i’  s ' V  V,' 
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[1  +  0.075  1  x  si8»  (v»bl)  (A— 6) 


Here  is  the  Reynolds  number  based  on  chamber  height  and  the 

corresponding  relative  flow  velocity  for  fixed  or  rotating  surfaces. 

The  factor  sign  (vRgL)  is  needed  to  give  the  forces  r.' ,  r"  ther  proper 

direction.  Thus,  we  append  the  factor  sign  (c*)  to  A'  and  the  factor 
$ 

sign  (wR,  -  c,  )  to  A",  both  here  and  in  the  first  order  calculations  of 
Appendix  2.  Dh  is  the  gland  hydraulic  diameter.  Eqs.  (A-3)  through 
(A-6)  can  be  solved  for  the  distribution  c,  of  azimuthal  velocities. 

In  particular,  the  asymptotic  velocity  (c*)  follows  from  (A-3)  when 

*  *  , 
ci  =  ci-i  1S  assuned. 


Appendix  2.  Coefficients  for  the  Perturbation  Equations 
Let  the  transverse  area  of  a  gland  (fig.  1)  be  f.  =  (h,  +  5;  )(i  in 
the  centered  position.  An  asterisk  on  any  variable  denotes  the 
undisturbed  (centered)  condition.  We  obtain  for  Eq.  (7)  the  following 
coefficients. 
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Here  t  is  the  ratio  of  specific  heats  and  R*  the  gas  constant. 

Appendix  3.  Results  for  a  Single-Cavity  Seal 
Sample  calculations  were  made  for  a  single-chamber 
straight-through  seal  with  6  =  0.25  mm.,  1  =  8  mm.,  h  =  3ram. ,  R,  =  15 
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cm.,  u  =  1000  rad/sec,  c  =  100  m/sec,  P0  =  1.5  atm,  Pa  =  1  atm.  For  a 
range  of  whirl'frequencies  from  0  =  0  to  (2  =  750  rad/sec,  the  resulting 
direct  and  quadrature  forces  are  shown  in  figure  A-l.  The  quadrature 
force  F, ,  which  is  the  one  of  importance  for  stability  considerations, 
is  seen  to  be  very  nearly  linear  with  0,  indicating  no  difference 
between  the  local  and  global  definitions  of  the  coefficients. 

There  is  is,  on  the  other  hand,  a  slight  curvature  in  the  Fd  line. 
The  values  calculated  for  Q  =  500  rad/sec  are  as  follows: 

K,x(N/m)  K*  x  (N/m)  Cn(N  sec/m)  Cyi  (Nsec/m) 

Local  (Eqs.  (15))  610  11350  +  13.7  -  0.95 

Global  (Eqs.  (16))  652  11480  13.98  -  1.02 

For  the  same  seal,  with  0  =  0  throughout,  figure  A2  shows  the 
effect  of  parametric  variations  of  the  friction  coefficients  (A '  and  A" 
varied  simultaneously)  at  various  shaft  rotation  speeds. 
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Tabic  2.  K  Calculated  vs.  Experimental  (u>  •  0) 

**  rot  it 

g 

CA  Exper.  Calcula. 

Run  f  Seal  Type  JC  Pa^Po  (m/s)  u  (x  10^  N/m) 


1 

S,  0,0,0 

17 

.66 

38.4 

.92 

.75 

.611 

2 

S, o', 0,0 

17 

.32 

52.9 

.92 

2.57 

2.091 

• 

3 

S, 1,0,0 

17 

.32 

68.3 

1.02- 

1.57 

1.45 

4 

S, 1,0,0 

17 

.66 

33.4 

1.02 

.27  * 

0.274 

5 

S, 1,0,0 

17 

.66 

48.5 

1.02 

.423 

0.411 

6 

S, 1,0,0 

17 

.56 

39.0 

1.02 

.457 

0.410 

7  * 

S, 1,0,0 

17 

.79 

38.2 

1.02 

.218 

0.224 

8 

S, 0,0,0 

17 

.49 

63.5 

.92 

1.89 

1.586 
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17 
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54.3 
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2.22 

1.747 

10 

S, 0,0,0 

17 

.79 

15.6 

.92 

.184 
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17 
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.92 
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1.605 

12 
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17 

.4* 
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.92 
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13 

S,1.0,0 

17 

.49 
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1.02 
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1.098 

14 

S.1,0,0 

17 

.49 
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1.02 

.57 
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. 

13 
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17 
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.92 

1.9 
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16 
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17 
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27.6 
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0.937 

17 
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18. 3Z 
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(m/a) 

(rad/s) 

u 
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23 
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Fig.  1  Geometry  for  labyrinth  seal  analysis 


Fig.  2  Direct  (F  )  and  quadrature  (F  )  forces 
P  q 

due  to  rotor  eccentricity 


Fig.  3  Two  different  definitions  of  spring  and 
damping  coefficients. 

(a)  Local,  (b)  Global 


CHAMBER  NUMBER  1 

Fig-  5  Calculated  vs.  experimental  axial  pressure  distribution 


CHAMBER  NUMBER  1 

Fig.  6  Calculated  vs.  experimental  circumferential  velocity 
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Fig.  7  Calculated  vs.  experimental  cross  spring 
coefficient  Kyx  (u>  ■  0) 
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Fig.  8  Entry  swirl  parameter  and  calculated  and 
experimental  cross  spring  coefficient  Kyx 


Fig.  A-l  Direct  and  Quadrature  Forces  in  a  Short  Seal 


Fig.  A-2  Effect  of  Friction  Factor  Variation  on  Cross-Stiffness 
in  a  Short  Seal.  Inlet. Swirl  Velocity  C  •  100  m/sec. 


